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1 Introduction 

Throughout the paper, we denote the real number field by M, the set of all m x n 
matrices over the quaternion algebra 

HI = {oo -I- aii + a2j -I- a^k \ P = = —1 , oq, oi, 02, 03 € R} 

by and by the set of all m x n matrices over H with a rank r. Let 

M (n, H) be the ring of n x n quaternion matrices. For A G the symbols 

A* stands for the conjugate transpose (Hermitian adjoint) matrix of A. The 
matrix A = (oy) G is Hermitian if A* = A. 

In the past, researches into the quaternion skew field had more a theoretical 
importance, but now a growing number of investigations give wide practical 
applications of quaternions. In particular through their attitude orientation, 
the quaternions arise in various fields such as quaternionic quantum theory [T] , 
fluid mechanics and particle dynamics ma, computer graphics [1], aircraft 
orientation [5], robotic systems [5], life science [ZlIB] and etc. 

Research on quaternion matrix equations and generalized inverses, which 
are usefulness tools used to solve matrix equations, has been actively ongoing 
for more recent years. We mention only some recent papers. Yuan, Wang and 
Duan [a derived solutions of the quaternion matrix equation AX = B and their 
applications in color image restoration. Wang and Yu m studied extreme ranks 
of real matrices in solution of the quaternion matrix equation AXB = C. Yuan, 
Liao and Lei El obtained the expressions of least squares Hermitian solution 
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with minimum norm of the quaternion matrix equation (AXB, CXD) = (E, F). 
Feng and Cheng m gave a clear description of the solution set to the quaternion 
matrix equation AX — XB = 0. Jiang and Wei m derived the explicit solution 
of the quaternion matrix equation X — AXB = C. Song, Chen and Wang m 
obtained the expressions of the explicit solutions of quaternion matrix equations 
XF—AX = BY and XF—AX = BY. Yuan and Wang [TS] gave the expressions 
of the least squares ? 7 -Hermitian solution with the least norm of the quaternion 
matrix equation AXB + CXD = E. Zhang, Wei, Lia and Zhao derived [12] the 
expressions of the minimal norm least squares solution, the pure imaginary least 
squares solution, and the real least squares solution for the quaternion matrix 
equation AX = B. 

The definitions of the generalized inverse matrices have been extended to 
quaternion matrices as follows. 

The Moore-Penrose inverse of A € denoted by Al, is the unique 

matrix X G satisfying the following equations 

1) AXA = A; 2) XAX = X; 3) (AX)* = AX; 4) (XA)* = XA. 

For A G with k = Ind A the smallest positive number such that rank A^+^ 

rankA^, the Drazin inverse of A, denoted by A^, is defined to be the unique 
matrix X that satisfying Eq. 2) and the equations 

5) AX = XA; 6) A^^+^X = A'^. 

In particular, when IndA. = 1, then X is called the group inverse of A and is 
denoted by X = A®. If IndA = 0, then A is nonsingular, and A^ = A’l’ = A^^. 

Cline and Greville m extended the Drazin inverse of square matrix to rect¬ 
angular matrix, that has been generalized to the quaternion algebra as follows. 

For A G and W G the W-weighted Drazin inverse of A with 

respect to W is the unique solution to the following equations 

7) (AW)'=+iXW = (AW)'=; 8) XWAWX = X; 9) AWX = XWA, 

where k = max{/nd(AW),/n(i(WA)}. 

The Drazin inverse and weighted Drazin inverse has several important appli¬ 
cations such as, applications in singular differential and difference equations m, 
signal processing m, Marckov chains and statistic problems uniin], descrip¬ 
tor continuous-time systems |22j . numerical analysis and Kronecker product 
systems [23], solving singular fuzzy linear system [24l]^, constrained linear 
systems [52] and etc. 

Cramer’s rule for the W-weighted Drazin inverse solutions, in particular, 
have been derived in m for singular linear equations and in m for a class 
of restricted matrix equations. Recently, within the framework of the theory 
of the column-row determinants Song [28] has first obtained a determinantal 
representation of the W-weighted Drazin inverse and Cramer’s rule of a class of 
restricted matrix equations over the quaternion algebra. But in obtaining, he 
has used auxiliary matrices other than that are given. In |29j , we have obtained 
new determinantal representations of the W-weighted Drazin inverse over the 
quaternion skew field without any auxiliary matrices. 

An important application of determinantal representations of generalized 
inverses is the Cramer rule for generalized inverse solutions of matrix equations. 
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But, when there is a need for a W-weighted Drazin inverse solution? Con¬ 
sider for example the following matrix equation, AiX = D. If Ai is rectangular 
and we can represent it as Ai = WAW, where WA and AW are quadratic 
and singular, then its W-weighted Drazin inverse solution is needed. 

In the paper we investigate analogs of Cramer’s rule for W-weighted Drazin 
inverse solutions of the following matrix equations over the quaternion skew 
field H, 


( 1 ) 

( 2 ) 

(3) 


WAWX = D, 
XWAW = D, 
W1AW1XW2BW2 = D. 


The paper is organized as follows. We start with some basic concepts and 
results from the theory of the row-column determinants in Subsection 2.1. In 
Subsection 2.2, we give determinantal representations of the Moore-Penrose and 
Drazin inverses for a quaternion matrix. Determinantal representations of the 
W-weighted Drazin inverse and its properties we consider in Subsection 2.3. 
In Subsection 3.1, we give the background of the problem of Cramer’s rule for 
the W-weighted Drazin inverse solution. In Subsection 3.2 we obtain explicit 
representation formulas of the W-weighted Drazin inverse solutions (analogs of 
Cramer’s rule) of the quaternion matrix equations (IT]131). In Section 4, we give 
numerical examples to illustrate the main result. 

2 Preliminaries 

2.1 Elements of the theory of the column and row deter¬ 
minants 

The theory of the row-column determinants over the quaternion skew field has 
been introduced in [50115^ , and later it has been applied to research generalized 
inverses and generalized inverse solutions of matrix equations. In particular, 
determinantal representations of the Moore-Penrose [331131] and explicit repre¬ 
sentation formulas for the minimum norm least squares solutions of some quater¬ 
nion matrix equations [35] , and determinantal representations of the Drazin [36j , 
and W-weighted Drazin inverses [29] have been obtained by the author. Song 
at al. derived determinantal representation of the generalized inverse g [37], 
Bott-Duffin inverse [55] and the Cramer rule for the solutions of restricted ma¬ 
trix equations [53], and the generalized Stein quaternion matrix equation m, 
etc. 

For A = (aij) £ M(n,]HI) we define n row determinants and n column 
determinants as follows. 

Suppose Sn is the symmetric group on the set n}. 

Definition 2.1 The i-th row determinant of A = (oij) £ M(n,]Hl) is defined 
for all i = l,n by putting 

rdetiA = (~1) + l ■ -^ 41 + 11 * ■ • ■ ®*Jsr»fcr + l ■ • ■ ®*»!r + ir*fcr ’ 

creSn 

^ + l ■ • ■ ikx+lf) ifk^^k^ + l • ■ • ik2+l‘z) ' ' ' (^fcr^fcr-t-1 ■ ' ■ ikT+lf) i 

with conditions iu^ < ik3 < ■ ■ ■ < ikr 0,'nd ikt < ikt+s for t = 2 , r and s = l,lf 
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Definition 2.2 The j-th column determinant of A = (aij) G M(n,IHI) is de¬ 
fined for all j = l,n by putting 

cdetj A = (~1) ^jkr.jkr+lr ' ' ' ®ifcr + l*fcr ' ' " ' ' ®ifci+lifci 

rdSn 

~ {jkr+lr ■ ■ ■ jkr + ljkr) ■ ■ ■ 0^2+^ ' ' ' jk2 + ljk2) (jfel+il ■ • ■ jfcl + ljfelj) I 


with conditions, jk 2 < jk^ < ■ ■ ■ < jk^ <knd jk^ < jkt+s for t = 2,r and s = IJt- 

Suppose denotes the submatrix of A obtained by deleting both the ith row 
and the jth column. Let a.^ be the j-th column and a^. be the *-th row of 
A. Suppose A j (b) denotes the matrix obtained from A by replacing its j-th 
column with the column b, and A^. (b) denotes the matrix obtained from A by 
replacing its i-th row with the row b. 

The following theorem has a key value in the theory of the column and row 
determinants. 

Theorem 2.1 1,W1| / If A = (aij) G M (n, H) is Hermitian, then rdeti A = • • • = 

rdet„A = cdetiA = • • • = cdet„A G R. 

Since all column and row determinants of a Hermitian matrix over H are equal, 
we can define the determinant of a Hermitian matrix A G M(n,]HI). By defi¬ 
nition, we put, detA := rdeti A = cdeti A, for all i = l,n. The determinant 
of a Hermitian matrix has properties similar to a usual determinant. They are 
completely explored in [sniisi] by its row and column determinants. They can 
be summarized by the following theorems. 

Theorem 2.2 If the i-th row of a Hermitian matrix A G M(n, H) is replaced 
with a left linear combination of its other rows, i.e. a^. = Ciai^, Ck&ik., 

where c; G HI for all I = l,k and {i, ii} C In, then 


rdeti Ai. (ciai^, Ck^i^.) — cdeti Ai, (ciai^, CkB-i^.) — 0. 

Theorem 2.3 If the j-th column of a Hermitian matrix A G M(n,]HI) 

is replaced with a right linear combination of its other columns, i.e. a j = 
a.jiCi -I- ... -I- a.jj,Cfc, where c; G H for all I = l,k and {j,ji} C J„, then 

cdetj A,j (a.jjCi a,j^Ck) = rdetj A.^ (a.jjCi a,j,_Ck) = 0. 

The determinant of a Hermitian matrix also has a property of expansion along 
arbitrary rows and columns using row and column determinants of submatrices. 

We have the following theorem on the determinantal representation of the 
inverse matrix over H. 


Theorem 2.4 1,W1| / The necessary and sufficient condition of invertibility of 

A G M(n,]HI) is ddetA ^ 0. Then there exists A~^ = (LA) ^ = (RA) 
where 


{LA)~^ = (A*A)~^ A* 


1 

ddetA 


t'hii 

Li2 

\Llra 


L 21 

L 22 

lb2n 


Lnl\ 

lb„2 

Hnn / 


(4) 
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{RA)~^ = A* (AA*)-^ 


1 

ddetA^ 


11 


•21 


(5) 


^12 


^22 


Rnl\ 

Rn2 


\Rln R2n ■■• Rnn/ 


and 

hij = cdetj(A*A),j (a*J , K ^ = rdetj(AA*)i, (a*) , 
for i,j = l,n, and ddetA = det(AA*) = det(A*A). 


2.2 Determinantal representations of the Moore-Penrose 
and Drazin inverses over the quaternion skew field 

We shall use the following notations. Let a := {ai,..., ak\ C m} and 

/3 := {/3i,..., /3fc} C {1,..., n} be subsets of the order 1 < fc < min {m, n}. By 
A^ denote the submatrix of A determined by the rows indexed by a and the 
columns indexed by /3. Then A“ denotes the principal submatrix determined by 
the rows and columns indexed by a. If A G M (n,]HI) is Hermitian, then by |A“| 
denote the corresponding principal minor of det A. For 1 < fc < n, the collection 
of strictly increasing sequences of k integers chosen from {1,... ,n} is denoted 
by Lk,n := {a : a = (ai,..., ak), 1 < ai < ... < a/c < n}. For fixed i G a and 
J G /?, let ^ ^ G Oi}, Jr, n f J G Lr,n^ j G /?}■ 

Denote by a* and a* the j-th column and the i-th row of A* and by a*'™^ 

and a|™^ the j-th column and the i-th row of A™, respectively. 

The following theorem give determinantal representations of the Moore- 
Penrose inverse over the quaternion skew field H. 

Theorem 2.5 If A G then the Moore-Penrose inverse A+ = 

(a^-) G possess the following determinantal representations: 


I 0£Jr,n{i} 

^^J = - 


E cdeD ((A*A) Eay)) ^ 


E (A*A) ^ 


3eJr. 


( 6 ) 


or 

rdetj ((AA*)j. (a* )) “ 

-I- _ «gF,mD} _ , , 

E l(AA*)“| ■ 

aG/r, m 

for all i = l,n, j = 1, m. 

Proposition 2.1 l2Uf If Ind{A) = k, then A^ = A^(A^^+^)+A^. 

Denote by a.g and at, the s-th column of (A^^+^)*A^ =: A = (a^) G 
and the t-th row of A^(A^^+^)* =: A = (a^) G respectively for all s,t = 

l,n. Using the determinantal representations of the Moore-Penrose inverse 
and and Proposition 12.11 the following determinantal representations of the 
Drazin inverse for an arbitrary square matrix over H have been obtained in [33] . 
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Theorem 2.6 \SSfl //A S M with IndA = k and rank = rankA*^ = 

r, then the Drazin inverse A^ possess the determinantal representations, 


E E cdett ((A^'^+i)* (A^'^+i) Ja.j)) 


JO * —1 /36J^r, n{t} 

^ij — 


x: (A2fc+i)*(A2fe+i)^ 

/3eJr.n 


and 


„D _ 

O-ij — 


ti Y. rdet«((A2'=+i(A2'=+i)*) (a,.))“)a 

s=l \ae/r,n{s} y .S / J 


(k) 

sj 


Y |(A2fe+i (A2fe+i)*) 

aG/r n 


( 8 ) 


(9) 


In the special case, when A S M (n, H) is Hermitian, we can obtain simpler 
determinantal representations of the Drazin inverse. 


Theorem 2.7 fddf //A S M (n, H) is Hermitian with IndA = k and rank 
rank A^ = r, then the Drazin inverse A^ = (a^) G possess the following 

determinantal representations, 


D 

— 


E cd*((AWi) .«))? 

3GJr,n{i} 


/3GJr 


E (A'^+i) i 


( 10 ) 


or 




E rdetj ((A'=+i)j. (afM 
E l(A'=+i)g| 

Ot ^ /r, n 


( 11 ) 


2.3 Determinantal representations of the W-weighted Drazin 
inverse 

Definition 2.3 For an arbitrary matrix over the quaternion skew field, A G 
jjmxn^ we denote by 

TZr{A) = {y G H™ : y = Ax, x G H”} , the column right space of A, 

J\fr{A) = {y G H” : Ax = 0}, the right null space of A, 

TZi{A) = {y G H" : y = xA, x G in’"}, the column left space of A, 

AfriA) = {y G H"* : xA = 0}, the left null space of A. 

We introduce some mathematical background from the theory of the W- 
weighted Drazin inverse [miiijiig that can be generalized to H. 

Lemma 2.1 Let A G andW G withk = max{/nd(AW), Jnd(WA)}. 
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Than we have: 


(a) Ad.w = A((WA)^))" = ((AW)^))"A; 

(b) Ad,wW = (AW)^; WAd,w = (WA)^; 

(c) Ad,w = {(AW)'= [(AW)2 '=+i]+ (AW)'=} W+; 

Ad,w = W+ |(WA)'= [(WA)2 '=+i]+ (WA)'=|; 

(d) WAWA£;,W = P-Rr((WA)'“),A/',.((WA)'“)) Ad,wWAW = P7Jj((aW)'“),M((AW)'=) i 

where P 7 j,^((wA)'').A^r((WA)'^) *5 the projector on 7?.r.((WA)^) along Mr{(WA.)'^), 
and P 7 ?,,((AW)''),A'i((AW)'=) ts the projector on TZi{{AW)^) along Afi{{AW)^). 

In particular, the point (a) of Lemma 12.11 due to Cline and Greville m is 
generalized [28] to H. Using this proposition, we have obtained [29] the following 
determinantal representations W-weighted Drazin inverse. 

Denote WA =: U = (u^) e and AW =: V = (vij) G 

Due to Theorem 12.61 we denote an entry of the Drazin inverse by 


D,1 

u,, = 


t 4'^ E cdeU ^ (u.,)) ^ 

t = l /3eJr,n{t} ^ ^ 


^ (U 2 fc+ 1 )* (U 2 '=+l) P 

3eJr,n 


or 


D,2 

= 


ti E rdeu((u2'=+i(U2'=+i)*) (u,.))“)u 

5=1 \ae7^,„{s} / .s / J 


(k) 

sj 


|(U2'=+l (U2fc+1)*) 

Ot^lr n 


( 12 ) 


(13) 


where u.^, and Ut. are the s-th column of (Xj2^+i)*U^ =: U = {uij) G and 

the t-th row of U*(U^*+^)* =: U = {iiij) G H”^”, respectively for all s,t = 1, n, 
r = rankU'=+i =rankU^ 

Then we have the following determinantal representations of A^i^w = G 


n 



d,W / D\(2) 

% =Z^«*9Ki) 

9=1 

(14) 

where 

71 



= E<p'4/ 

p^i 

(15) 

for all 1, / = 1, 2, and 
Similarly using V 

u^j’^ from ([T^ and from (fT^. 

= {vij) G 



m 

(16) 


9=1 
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where the first factor is one of the four possible equations 


p=i 

for all I, f = 1,2, and an entry of the Drazin inverse is denoting by 


E 




D 1 t-1 /3GJr,m{t} 

= - 


E cdett ((V2'=+1)* (V2'=+') , (V,)) 


j2 ^ 

^Jr,m 


or 


D,2 

= 


t( E rdet«((v2^+i(V^''+')1 (v..))sl 

S=1 \ae/,-,™{s} / .s / J 


,ik) 

sj 


I (V2'=+l)*) 

aG/r m 


(17) 


(18) 


(19) 


where v.s and vj, are the s-th column of (V^*+^)*V^ =: V = (vij) € 
and the t-th row of V^(V^^+^)* =: V = {vij) G respectively for all 

s,t = l,Tn, r = rankV^+^ = rankV*. 

The point (c) of Lemma 12.11 due to [53] has been generalized to H in [33] . 
Using this proposition, we have obtained the following two determinantal rep¬ 
resentations of the W-weighted Drazin inverse. 

Theorem 2.8 129f Let A G ]]j™xn g with k = Ind{AW) and 

r = rank(AW)^+^ = rank(AW)*. Then the W-weighted Drazin inverse of A 
with respect to W possesses the determinantal representations 


d,W 

a^J = 


m 


X; E rdett 

aGlr-^ m{i} 


(^(Y2k+1 ^Y'^k+iy'^ 

IJv.,)) 

a E I'detj 1 

aG/rl,,^{i} 

((WW*)^. (wO) 

E 

y/-2k-\-l ('Y2fe+1)*^ a 

1 E l(ww*)“| 


OL^ir,m OL^Iri,n 


Ot 

CX. 


( 20 ) 


and 


d,W 

% = 


E E 

i=l/3GJri,m{*} 


cdeti ((W*W) J (w.t)) ^ X) cdett 

Jr, Ti{i} 


^(uj) I 




(w*w)^ ^ ((U2fe+i)*u2fe+i) 


3GJr, 


(21) 


where V = v'=(V2'=+i)*, W = V'^W*, and U = (U2'=+l)*U^ W = W*U'=. 

In the special cases, when AW G jjmxm WA G are Hermitian, 

we can obtain simpler determinantal representations of the W-weighted Drazin 


inverse. 













Theorem 2.9 129^ If A G W G and AW G is Hermi- 

tian with k = max{Jnd(AW),/n(i(WA)} and rank(AW)^+^ = rank(AW)^ = 

r, then the W-weighted Drazin inverse Ad^w = ^ jjmxn respeet to 

W possess the following determinantal representations: 


d,W 

a, o’ = 


/3eJ, 


E^^cdet.((AW)^f (v.,)) 


E (AW) 

ySeJr, m ' 


k+2 p 


( 22 ) 


Theorem 2.10 If A G H™""", W e and WA G H"''" is Hermi- 

tian with k = inax{Ind{AW), IndfWA)} and rank(WA)^+^ = rank(WA)^ = 


where v ^- is the jth column of'V = (AW)^A for all j = l,m. 


r, then the W-weighted Drazin inverse Adpv = 

W possess the following determinantal representations: 

E rdet,((WA)^H2(-WA 

d,W _ aeJr.nb'} _ 


E (WA) 

Oi^lr, n 


k+2 


with respect to 


(23) 


where u^. is the ith row ofXJ = A(WA)^ for all i = l,n. 


3 Cramer’s rule for the W-weighted Drazin in¬ 
verse solution 


3.1 Background of the problem 


In ^7] Wei has established a Cramer’s rule for solving of a general restricted 
equation 

WAWx = b, xG7^[(AW)'=l] , (24) 

where A G W G with Jnd(AW) = ki, Ind (WA) = fca and 

rank(AW)^i = ri, rank(WA)^2 = r 2 . He proofed that the restricted matrix 
equation (IMl) has a unique solution, x = Ad^w^, and presented its Cramer’s 
rule as follows, 


Xj = det 


fWAW{j —s- b) 

V Vi(j^O) 




fWAW 

V Vi 



(25) 


where Ui G G are matrices whose columns form bases 

for M(^(WA)+ and A/'((AW)'=*), respectively. 

Recently, within the framework of a theory of the column and row determi¬ 
nants Song has considered the characterization of the W-weighted Drazin 
inverse over the quaternion skew and presented a Cramer’s rule of the restricted 
matrix equation. 


W 1 AW 1 XW 2 BW 2 = D, (26) 

7^dX) c TZr {{AW+^) , A;(X) d Mr ((W 2 B)'=’^) , (27) 

7^^(X) c III ((BW 2 )'==) , M(X) D Ml ((WiA)'=i) , 
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where A S Wi S B e W 2 S and D e with 

ki = max{/nd(AWi),/nd(WiA)}, fca = max{/nd(BW 2 ),/nd(W 2 B)}, and 
rank(AWi)^i = si, rank(BW2)*^^ = S2- 
He proofed that if 

D e 7^, ((WiA)'=i, (W2B)'=^) , D e 7^^ ((AWi)'=i, (BW2)'==) 

and there exist auxiliary matrices of full column rank, Li S M* S 

L 2 G M 2 G with additional terms of their ranges 

and null spaces, then the restricted matrix equation (PH|) has a unique solution, 

X = Ad,WiDBd,w2- 

Using auxiliary matrices, Li, Mi, L 2 , M 2 , Song presented its Cramer’s rule by 
analogy to (1^ . 

In this paper we have avoided such approach and have obtained explicit 
formulas for determinantal representations of the W-weighted Drazin inverse 
solutions of matrix equations by using only given matrices. 


3.2 A Cramer’s rule for the W-weighted Drazin inverse 
solutions of some matrix equations 

Consider the following restricted matrix equation, 

WAWX D, (28) 

7^dX) c TZr {{AWf) , Mi{X) D M ((WA)'=) , (29) 

where A G W G with k = max {Jnd(AW),/nd (WA)}, and 

D G mu'^p. 

Theorem 3.1 //D C TZr ((AW)^) and D D A/j ((WA)^), then the restricted 
matrix equation IgdI) has a unique solution, 

X = Ad,ivD, ( 30 ) 

which possess the following determinantal representations for all i = l,m, j = 

i) 


Xij — 


X; E cdeU ((W*W) j (w.t)) ^ E cdett 

i = l/3GJri,m{i} ^aeJr, n{t} 


U2fe+l)*U2fe+l^ (dj-)) 




E (w*w)^ E ((U2fe+i)*u2fe+i) 


3eJr, 


(31) 


where d j is the j-th column 0/ D = UD = (U2'=+i)*U'=D, U = WA, W = 
W*U^ and r = rank(WA)'=+i = rank(WA)'=. 

a) 

m 

( 32 ) 

9=1 
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where can he obtained by and AD = R = (vqj) G 

in) If AW G is Hermitian, then 

E cdet.((AW)'=f (f.)) ^ 

Xij — 

E 

P&Jr-, 

where f.^ is the j-th column o/F = VD = (AW)^AD. 

Proof. The proof contains two parts. We first shall establish that the unique 
solution of ((^ can be represented as (111(1 . By the definition of the right range, 
we have 

D = (WA)'=Y 

for some matrix Y G IHI"^^. It follows that 

7^^(D) C TZr {{WAf) . 

Then by Lemma [2.II (dl. 

WAWAd^vvD = P7?,,.((WA)''),A/'r((WA)'')D = D. 

It means that dsni is a solution of (EHD and satisfies the restricted conditions 

dMl). 

Now we prove the uniqueness of dsni). Let Xo is a solution of (l28l) . Then it 
satisfies the restricted conditions (EHl) . and 

= Ad^rvWAWXo = P7?,,.((WA)‘^),A/'r((WA)'“)Xo = Xq. 

To derive a Cramer’s rule (I1TT|1 . we use the determinantal representation GB) 
for Ad^w- Then 


(AW) 


fc+2 /3 




E 

S = 1 


E d,W , 

flzs dsj = 

S^l 

t E cdeL(W*W) E cdeL((U2'=+i)*U2'=+i) (u., 

t = l /3eJr, . mid 0eJr.n{t} ^ ^ * 


E (w*w)^ E ((U2fe+i)*u2'=+i) ^ 

ySeJri.m /SGJr.n. 


(34) 


Denote D = UD = 


(U2fc+i)*u'=D, where D = (dsj^ G Since 


P 

E! ^ sdsj = d.j, 

S=1 

where d.j is the j-th column of D, then (1311) follows from (IMl) . 

Similarly, we derive the analogs of Cramer’s rule (l32ll and (IMl) by using the 
determinantal representations for the W-weighted Drazin inverse ((T5)) , and (1^^ , 
respectively. □ 
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Remark 3.1 In the complex case, i.e. A G W G and D G 

C^'^P, we substitute usual determinants for all corresponding row and column 
determinants in and hSS\] . 

Note that in the case in), the condition AW G be Hermitian is not 

necessary, then in the complex case iSS\) will have the form 

E |((AW)^f (f.,)) ^ 

Xij — 

E 

P&Jr-, 

where i,j is the j-th column o/F = VD = (AW)^AD. 

Now, consider the following restricted matrix equation, 

XWAW = D, (35) 

TZiiX) c III {{AWf) , A;(X) d Nr ((WA)'=) , (36) 

where A G W G with k = max {Jnd(AW),/nd (WA)}, and 

D G 

Theorem 3.2 IfDcTZi ((AW)^) and D D N ((WA)^), then the restricted 
matrix equation has a unique solution, 

X = DAd,vy, (37) 

which possess the following determinantal representations for i = l,q, j = 1, n, 

i) 

Xij = 

E E rdeti((v2'=+i(V2'=+i)*) (d,.))“ E rdet, ((WW*) . (wE) “ 

E |(v 2 fe+i(V 2 fc+i)*) g| E l(ww*)g| 

Ot€lr,m Oi£lri,n 

(38) 


(AW)''+^ ^ 


where di, is the i-th row o/D = DV = DV'=(V2'=+i)% V = AW and r = 
rank(AW)'=+i = rank(AW)'=. 
ii) 

n 

Xi, =^kt{u^j)^'^\ (39) 

where (u^)^^) can be obtained by if73l) and DA = L = (In) G 
Hi) If AW G is Hermitian, then 


E rdet^ 

«6Jr,n{i} 

{{WA)';+\gr)) “ 

E 

OiGlr, n 

(WA)'=+" “ 



(40) 


where gi, is the i-th row of G = DA(WA)^ for all i = l,n. 
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Proof. The proof is similar to the proof of Theorem 13.11 □ 


Remark 3.2 In the complex case, i.e. A € W € and D £ 

C"^^, we substitute usual determinants for all corresponding row and column 


determinants in i38\) . iStA) . and {4C^ . Herein the condition WA € C"^" be 
Hermitian is not necessary, then in the complex case can be represented as 
follows. 


Xij — 


E 

a&Ir,n{j} 


|((WAi"+2(g,.)) 


E 

OL^Ir 


(WA) 


k+2 


where gi, is the i-th row of G = DA(WA)^ for all i = l,n. 


Now we consider the matrix equation ([^ with the constraints (EU. Denote 
ADB =: D = (dif^ € and VDU =: D = (dif) G where 

V := (AWO'^iA , U := B(W2B)'==. 


Theorem 3.3 Suppose D G A G Wi G with ki = 

max{/nd(AWi),/nd(WiA)}, where rank(AWi)'=i = si, and B G 
W 2 G with k 2 = max{/nd(BW 2 ),/nd(W 2 B)}, rank(BW 2 )'== = S 2 . 

//D G 7^'((WlA)'=^(W2B)'==), D G 7^^ ((AWi)'=N(BW 2 )'==), then the re¬ 
stricted matrix equation has a unique solution. 


X = A. 


d,WiDBd,w2i 


(41) 


which possess the following determinantal representations for all i = l,m, j = 
Tfq. 

i) 

m q 

(42) 

1=1 /=! 

where (v^) = is the Drazin inverse of V = AWi and can be 

obtained by and (u^j) = is the Drazin inverse 0 / U = W 2 B and 

can be obtained by m- 

a) If AW I G and W 2 B G are Hermitian, then 


E cdetJ(AWi)':r^(dB)K 


E (AWi) 

Jsi, m 


fci+2 /3 

p 


E (W 2 B) 
Ol^ I So , Q 


A;2~I”2 


(43) 


where 


E rdet, ((W 2 B)^"^+ 2 (dA)) 
cx^Is2 ,q{j} 


E (AWi) 

PeJsi.m 


fci+2 /3 
P 


E (W 2 B) 

OtGlso , Q 


^2 + 2 


(44) 


d",= 


^ rdet, ((W2B)^^=+2(dE) J G 


i = 1, n 


(45) 
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are the column vector and the row vector, respectively, d^. and d,j are the i-th 
row and the j-th column o/D for all i = l,n, j = l,p. 

Proof. The existence and uniqueness of the solution m can be proved similar 
as in ( [5S], Theorem 5.2). 

To establish a Cramer’s rule of (l26ll we note that we shall not use the deter- 
minantal representations (l2^ and (l22|) for dUD because corresponding determi- 
nantal representations of it’s solution will be too cumbersome. 

To derive a Cramer’s rule (1421) we use the sentence (a) from Lemma o 
Then we obtain 


X = ((AWi)^))" ADB ((WzB)^))" . (47) 

Denote ADB =: D = (dif'^ G V := AWi, and U := W 2 B. Then the 

equation (HZD will be written component-wise as follows 


p n 


d,Wl\i /ud,W2\ 


s=l 


p n / m 


dts 

s=l t=l \/=l / \/=! 


By changing the order of summation, from here it follows (H^ . 

ii) If A G B G WP^<i and AWi G and W 2 B G are 

Hermitian, then by Theorems 12.91 and 12.101 the W-weighted Drazin inverses 
Ad.Wi = ^ and B^^vVa = ^ posses the following 

determinantal representations respectively, 


E cdeCCAWi) ) (v.j;, ^ 

d,Wi _ 


E (AWi) 

3eJr, m 


fel+2 0 


(48) 


where v ^- is the j-th column of V = (AWi )^i A for all j = 1, m; 


E rdet, ((W2B))=+2(-^ a 


E 



^ ^ 2 , q 




(49) 


where u^, is the i-th row of U = B(W 2 B )^2 for all i = l,p. 
By component-wise writing dnj we obtain. 


Xij 




.d,W2 

^sj 


s=l \t^l 


(50) 
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Denote by the s-th column of VD = (AWi)'=i AD =: D = (4) e 
for all s = l,p. It follows from = d.g that 

t 


E d,Wli 

dts = 2^ 


E Cdet, ((AWi)^f Ev.t)) 


m{^} 


t=l 




E (AWi)'=^+^ ^ 


' dts — 


E Ecdet,((AWi)Y'(v, 0 ) 


l3eJsi,, 


E (AWi)"^+" ^ 


E cdeD((AWi)Y"(d..)) 


E (AWi)''^+^ ^ 

pGJsi 


(51) 


Suppose Gg. and e.g are respectively the unit row-vector and the unit column- 
vector whose components are 0, except the s-th components, which are 1. Sub¬ 
stituting m and dSl) in Q, we obtain 


Xij — 

E cdet.((AWi)'=(+2(d.g)) ^ E rdet, ((W 2 B)^"=+^(ug.,^ „ 
PeJs^,m{^} ^ V // ae/.2,,D} ^ 


E 


E (AWi)''^+^ ^ 


E (W2B) 

a&Iso. 


k2+2 


Since 


d. g — ^ ( G tdtsj ^s. — ^ ^ j ^ ^ digUsi — dti , 


(52) 


Z=1 


then we have 


^^3 
p n q 




E E E E cdet. ((AWi)^f ^ (e.O 1 ^dts^g, E rdet, ((W 2 B);=+^(ez, 




r 


E (AWj 


\ki +2 /3 


E (W2B) 

aGiso, Q 


/C2 +2 


n q 




EE E cdetE(AWi)^fEe.i)) E rdet, ((W 2 B);=+^(e, 


“6^32,qiii 


\k2+2( 


E (AWi) 

fieJsi.m 


ki +2 ^ 

/3 


E (W2B) 

Q! e I So , a 


^2 +2 


(53) 
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Denote by 


n 

Y, cdet.((AWi)Y'(d,z)) ^ = ^ ^ cdet.((AWi)Y'(e,0) 

0eJsi.m{i} ‘ = 1 /3eJsi,m{»} 

the Z-th component of a row-vector = {d^,..., d{^) for all I = l,q. Substitut¬ 
ing it in (l5!?ll . we have 


E rdet,((W2B)^^^+2(eO) 


Z=1 


OiGls^ 


Xij — 


Aj} 


E (AWi) 


ki+2 f3 

P 


E (W2B) 

CK^Isn , 


/C2 + 2 


Since E then it follows (l4^ . 

1=1 

If we denote by 


dfr-= 


Ydti Y rdet, ((W2B)^^^+2(e,)j^ = ^ rdet, ((W2B)^^^+^(dt 

z=i ae/« 2 ,,{i} aei, 2 ,q{j} 


(54) 


the t-th component of a column-vector d® = {dfj, for alH = l,n and 

substituting it in (l5^ . we obtain 


E E cdetJ(AWi)V'(e.t 

* = l/3gJ.i,m{d ^ _ 

E (A*A)^ E l(BB*)^ 

/3gi/r^,Ti ^^3r2,p 


/3 jB 
/3 “tj 


Since E G.idjj = d®, then it follows (@3]). □ 


4 Examples 

In this section, we give examples to illustrate our results. 
I. Let us consider the matrix equation 


WAWX = D 

with the restricted conditions (1^ . where 


/O 

k 

1 

VI 


i 

1 

0 

-k 




0 i 
k 0 
1 0 



(55) 
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Then 


V = AW 


^ —k —j 0 z ^ 
-1-j i + k j l+j 
k 0 i 0 
\—i + k 1—j i i — kj 


U = WA : 


j ON 
0 fc 0 
.0 0 0 ; 


and rankW = 3, rankV = 3, rankV^ = rankV^ = 2, rankU^ = rankU = 2. 
So, IndY = 2, Ind\J = 1, and k = max{Jnd(AW), Ind{WA)} = 2. 

We shall find the W-weighted Drazin inverse solution of by it’s deter- 
minantal representation (1311) . We have 

^-1 i + k 0\ A 2 + 3j 0^ 

= I 0 -1 1 = 0 k 

000 / \o 0 0 ; 


0 0\ /I -2i - 3fc 0^ 

T5^* Tt5 _ 2 


(U'^)*=(2-3j -k ,(U'^) U®= 2i + 3fc 

0 0 / \ 0 
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0 


/ i-j-k 

-j \ 

f-k 

0 

j 

-k 

o\ 

1 

1 + 3i + 6j - 2fc 4i - 2fc , W* = 

—i 

0 

0 

\ ^ 

0 / 


1 

k) 


w*w = 


Since by (1^ 


(^ 

i 

-j 

7 ^ 


(-k 

l-2j 


—i 

2 

0 

-2k 

,w = W*V‘^ = 

0 

i + k 

0 

j 

0 

1 

0 

i 

1 + J 

0 

\-J 

2k 

0 

2 ^ 



-1 

oy 


a^ii = 
3 


E E cdeti ((W*W) Ew,)) ^ E cdet*(((U^)*U5) (d.i)) 

‘=i/3e/3,4{i} /3eJ2,3{t} ^ ^ 


E (w*w)^ E ((U 5 )*U 5 )^ 

I3&J2,3 


where 


^ cdeti ((W*W), (w.i)) I = 

/36/3,4{1} 

A i —j\ /k 

0 2 0 

.i 0 1 


i 

■i \ 

r 

-j 

/\ 

2 

-2k 

+ cdeti i 

1 

0=0 

2k 

1 j 

\o 

0 

2 / 


^ cdeti ((W*W) , (w. 2 )) ^ = -2j, 

/5e/3,4{i} 

^ cdeti ((W*W) Ew. 3 )) ^ = 0, ^ |(W*W) 

/3e/3,4{l} 


= 2 , 
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and 


^ cdeti(((U5)*U5) ^(d.i)) 

/3G J 2 , 3{1} 


, , ( i—j — k —2i — 3k 

cdeti ( i_^3. _^6j_2fc 14 


+ cdeti 


i—j—k 0 
0 0 


= -22 -J - k, 


E cdet2(((U5)*U^) 

/3eJ2.3{2} 

^ cdet3(((U^)*U5) (d,i))^=0, ^ \{{Vr^'' 

/3eJ2,3{3} ■ /3eJ2,3 




= 1 , 


then 


0 ■ (-22 -j-k) + (-2j) • j + 0 ■ 0 
xn =--= 1, 

0 • (-2 + 2j) + (-2j) • i + 0 • 0 


X 12 =- 


2 ■ 1 


= k, 


2j-{-2i-j-k) + {10i-4:k)-j + 0-0 , , . , 

X 21 = -- = l+ i + 7k, 

2j-(-2 + 2j) + il0i-4:k)-i + 0-0 „ 

^22 =--= -7 - 4j, 

10i-(-2t-j-fc)+j-j + 0-0 

2^31 =-— - - = 9.5 + 5j — ok, 


X 32 =- 

We finally get, 


2 ■ 1 

lOi- (-2 + 2j)+j-i + 0-0 


2-1 


1 


= —lOi + 9.5fc, 

,v 


k 

X = ( l + i + 7k -7-4j 
\9.5 + 5j-5fc -102 + 9.5fc/ 

2. Let now we consider the matrix equation 


W 1 AW 1 XW 2 BW 2 = D, 

with the constraints (EH), where 


(56) 


(57) 


k 0 i 0 

A = I -j fc 0 1 I , Wi = 

0 1 0 -fc, 


B 


Since the following matrices are Hermitian 
'-2 i ON 


(k -j 0 \ 
0 fc 1 
i 0 0 

\0 1 -kj 


,W2=|j 0 

0 1 



/ 2 

-i\ 


fc 

0 

'vJ 0 l^’ 

0 

j 


El 

0 / 


'0 -i 


V = AWi =\-i -1 0 , U = W 2 B = 2-1 0 


0 0 0 ; 


.2 0 - 1 ; 
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then we can find the W-weighted Drazin inverse solution of (1571) by it’s deter- 
minantal representation (1431) . 

We have 


ki = max{/n(i(AWi),/nd(WiA)} = 1, 
k 2 = max{/nd(BW 2 ),/nd(W 2 B)} = 1, 

and Si = rank (AWi) = 2, S 2 = rank (W 2 B) = 2. Since 


then 


We have 


;SeJ2,3 


/-13 

8i 

0\ 

/O 

— 3z —3i\ 

-8i 

-5 

0 

, (W2B)3 = 3z 

-3 0 

V 0 

0 

0/ 

\3i 

0 3 / 

(AWi) 

3 /3 
/3 

= 1 

E (W2Bf“ 

= -27. 


3 


2i+j -7 + k -5 +2k'' 

D = AWi ADBW 2 B = \ -I + k -5i-j -4i -2j I , 

0 0 0 


By (gS]), we can get 


36z — 9j 


-27 


9-9fc' 


Since 


= -27- 9fc , d^2 = -18z , d^3 = 9z + 3j . 


36i — 9j 8i 0^ 
(AWi)"i (d^i) = I -27-9A: -5 0 1 , 


0 4; 


then finally we obtain 


X; cdeti [(AWi )^3 (d^''^ 


xn = 


/3eJ2,3{i} 


E (AWi)"^ E (W2Bf 

0^J2,3 Ck€/2,3 


36i — 27 j —Ai + 3j 

-27 ^ 3 ’ 


Similarly, 


cdeti 


a;i2 = 


-27 8i 
-18f -51 1 


cdeti 


= 3 , Xi3 = 


9 — 9k 8i 
9i — 3j —5j —9 — 7k 


cdet2 


2^21 = 


-27 

-13 36i - 9j 

-8i -27-9k I -7-5k 


-27 9 

"-13 -27 

—8i —18i 


cdet2 


cdet 


-27 3 

-13 -9-9k''' 


X 22 = 


-27 


-2i 


2^23 = 


-8i 9i + 3j J 15 i-nj 


2^31 = 2^32 = 2:33 = 0. 


-27 
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So, the W-weighted Drazin inverse solution of dsa are 


, /-12i + 9j 3 -9-7k\ 

X = - -21 - 15fc -6i 15i - llj . 

^ \ 0 0 0 / 

Note that we used Maple with the package CLIFFORD in the calculations. 
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